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SECTION  1.  INTRODUCTION  AND  Sim  ARY 


i.l  INTRODUCTION 

In  protection  against  as  EMP  [1]  (electromagnetic  pulse),  ferromagnetic 
metallic  shielding  is  one  of  the  simplest  and  most  commonly  used  effective 
schemes  [2] . The  cf fectiveness  of  the  lading,  measured  by  the  ratio  of  the 
field  penetrated  across  the  shield  to  the  incident  field,  depends  of  course 
qq  the  EM  properties  and  the  geometry  of  the  shielding  material.  Generally 
speaking,  the  dominant  material  properties  are  its  conductivity  s and 
permeability  y,  and,  if  the  shielding  enclosure  has  its  radius  of  curvature 
much  larger  than  the  «ave  length  of  the  incident  field  and  is  free  of 
seams  and  cracks,  the  only  important  geometrical  factor  is  the  thickness 
of  the  shielding  plate. 

In  this  report,  we  investigated  the  shielding  problem  for  the  incident 
EMP  in  the  form  of  a cylindrical  TEM  wave,  such  as  a wire  carrying  a 
surge  arrestor  current  terminated  at  an  iron  or  steel  shielding  plane  [3] . 
Because  of  the  large  strength  of  the  incident  field,  the  nonlinear 
ferromagnetic  saturation  of  the  plate  plays  an  important  role  in  deter- 
mining the  peak  and  the  shape  of  the  transmitted  field.  This  problem, 
based  on  and  together  with  its  one-dimensional  plane-wave-incidence 
version,  is  solved  analytically  for  a constant  u case.  The  analytical 
results  are  then  used  to  partly  predict  and  to  interpret  the  numerical 
results  for  the  one-dimension  nonlinear  case,  obtained  by  using  a finite 
difference  code  DIFUSN,  and  to  help  predict  the  behaviors  for  the 
cylindrical  nonlinear  case. 

In  the  following.  Section  1.2  briefly  summarizes  the  results  found  in 
this  report;  Section  2 solves,  analytically  and  numerically,  the  one- 
dimension  plane-wave  problem;  based  on  this.  Section  3 solves  the  cylin- 
drical incidence  probietn. 


I 


As  to  the  system  of  units  la  this  resort,  rationalised  !$KS  Is  used. 

1.2  SUMMARY 

^wnwaMii  mfcn 

W'e  briefly  sussaarize  here  the  results  obtained  la  this  report.  Details  of 
them  are  given  in  the  subsequent  text. 

1.  For  wavelengths  large  compared  to  the  shielding's  radius  of 

curvature,  high-p  conducting  plates  shield  EM?  very  effectively. 

For  short  *iulsa  and  thick  slab,  the  diff used-through  transmitted 

2 3-1 

field  varies  as  *{<*>  y a d ) . This  gives  a peak  transmitted 

-13  ° r 

field  v3  x 10  that  of  the  incident  peak  value  for  a typical 

ui  * 3 x 10^  radian/ sec,  a "vlO"*,  rv  10'  mho/meter,  d'vkna.  The 

° r 2 
else  width  of  the  transmitted  field  varies  as  <*uod  , giving 

‘'-3.3  x 10  1 sec  for  the  typical  example.  This  makes  the  shielding 

better  for  the  higher  frequencies,  and  thus  substantially  shifts 

downward  the  transmitted  wave’s  frequency  contents. 

2.  The  non-linear  saturation  of  the  ferromagnetic  permeability, 
which  saturates  d3/dH  to  smaller  values  for  stronger  field, 
slightly  reduces  the  transmitted  field  but  leaves  virtually 
intact  its  time  shape.  This  is  caused  by  the  fact  that  for  the 
relatively  narrow  SMP  in  a relatively  thick  slab,  the  saturation 
disperses  and  distributes  more  evenly  the  diffused  field  and 
mitigates  the  build  up  of  its  local  peak  value,  in  surprising 
contrast  to  the  simple  intuition,  naively  extrapolated  from  the 
const3nt-u  case,  that  a smaller  permeability  admits  more  field 
in  a shorter  time,  an  extrapolation  valid  only  when  a wide  pulse 
saturates  the  whole  thin  slab. 

3.  The  cylindrical  problem,  with  a wire  carrying  a large  current, 
such  as  a surge  current  arrestor,  terminated  at  a shielding 
wall,  is  solved  approximately.  Its  results  are  simply  related 
to  those  of  the  one-dimensional  problem,  enabling  us  to  make 


use  ©f  cha  o»e~diast lsios&I  results  for  the  cylindrical  pr«r* 
diet  isos.  Stoughly  speaking,  the  mala  part  oi  the  cylindrical 
diffused  fields  are  a/p  ti®es  that  of  the  oae~d isaejis iona  1 
anm  for  p^a  and  go  to  sera  far  a-*©,  with  minor  deviations 
being  complicated  functions  of  the  various  parameters  of 
the  p rob lea. 

The  shays  -emits  are  found  either  analytically  or  numerically,  or  both. 
Fro®  their ■ we  can  safely  conclude  that  for  practical  SMP  shielding®  using 
the  highly  effective  high-y  conducting  plates,  the  ferromagnetic  saturation 
only  slightly  enhances  the  shielding  effectiveness.  A side  result  is  that 
the  presences  of  holes,  crack®,  or  seams  at  the  shielding  plate  probably 
constitute  sore  important  nodes  of  penetration  for  the  incident  EM?. 
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SECTIOK 


A,  * 


THE  OSE-DIKESSION  SLAB  PROBLEM 


2.1  ANALYSIS  FOE  THE  CONSTANT  u CASE 


Consider  the  one-diaer.sional  problem  depicted  by  Figure  1.  A plane  electro- 
magnetic wave  with  fields 


Binc(t.i) 


iJaC(t,s)  - H0f(t  - |) 


. _ „inc,.  . _ -4/0 

1 (t’2)  * sx  (t>2)  * Bo!r 


■ Cc 


f la] 


(lb) 


is  incident  from  z * ■*»  normally  upon  a slab  of  thickness  d ac  position 
2 » D.  The  medium  to  the  left  and  to  the  right  of  the  slab  is  uniform 

and  has  a dielectric  constant  t , a permeability  u , and  a velocity  of 

-1/2  0 0 
light  e = (y  c ) {i).  The  slab  itself  is  also  uniform  and  has  a 

v OO 

dielectric  constant  e **  e e , a permeability  y ® y u , and  a conductivity 

TO  TO 

0 such  th3t 


0 >>  ye,  ye 

’ o 


(2) 


or  equivalently 


3 

If 


ji 1 

3t ' 


where  w is  the  angular  frequency  of  the  frequency  of  in*-ore«sr  \’n rice 

that  this  condition  of  high  slab  conductivity  is  assumed  throughout  this 

7 ~il 

report.  In  typical  cases,  we  have  3 v 10  mho/m  for  steel,  s v 10  Farad/m 
and  w E 10  radian/sec  for  the  incident  EMP;  thus  (2)  or  (2')  is  amply 
satisfied.  The  problem  is  to  find  the  electromagnetic  fields  everywhere, 
especially  in  and  transmitted  through  the  slab. 


Figure  1.  The  One-Dimensional  Problem 


2.1.1  General  Solution 


Being  a one-dimensicn  problem,  the  only  fields  are  plane  waves  with  a 
magnetic  field  in  the  y-direction  and  an  electric  field  in  the  x-direction. 
In  the  region  z < G,  the  reflected  fields,  in  addition  to  the  incident 
ones  given  by  (1),  are 


Hrfif(t,z)  = Ho  g(t  + f)  (3a) 

Eref(t,z)  * g(t  + (3b) 

In  z > 0,  the  transmitted  fields  are 

Htrans(tjZ)  = H T(t  _ iszllj  (4a) 

o c 

BCran3(t,z)  - T(t  - ^~L)  (4b) 

0 


finally.  In  the  slab  0 < i < d the  high  conductivity  condition  (2) 
combines  with  the  Maxwell  Equations  to  give  a simple  diffusion  equation 

3 2 , 

— r H(t » z)  - up  a — li(t,z)  * 0,  0 < z < d (5) 

3 £-  r u o t 

z 

for  the  magnetic  field  H(t,z)  to  which  the  electric  field  is  related  by 

E(t,z)  - ~~  |-  H(t, z)  (6) 

O oZ 

Now  to  find  the  fields,  we  merely  have  to  solve  (5)  and  (6),  subject 
to  the  boundary  conditions  that  require  continuous  magnetic  and  electric 
fields  at  both  slab  surfaces  z*0tdz“d.  In  terms  of  the  magnetic 
field  in  the  slab,  these  conditions  are  (Appendix  A) 


6 


4* 


1 , 

it 


li 


“ H<t,0)  oH(t,0)  - -2-JS 


0)  » -2 ^j~  R oa  f (t) 
~o 


H(t ,d)  + yS  0H(t,d>  - 0 


(7a) 


(7b) 


Thus,  Che  problem  reduces  to  solving  (5)  and  (7). 

To  express  the  field  in  a convenient  and  simple  form,  we  make  use  of  the 
Laplace  transform 


it 


H(s,z)  ■ f H(t,z)  e st  dt  = L[H(ttz)]  (8) 

JQ 


Then  from  (5)  and  (7),  the  resulting  transformed  fields  in  0 < z d 
are  (Appendix  B) 


H(s,z)  « 2Ho  f (s)  [A(s>  2 + B(s)  e*-**®*  z]  (9a) 

E(s,z)  * -2Hq  f(s)Vf  [-A(s)  e-^2  + B(s)  z]  (9b) 


where 


(9c) 


In  particular,  this  gives  the  transmitted  magnetic  field 


/ 


and  the  total  magnetic  field  at  the  incidence  side  of  the  slab 


Hq[!(s)  + g 


3)]  - H(s,0) 


2Hq  f (s) 

(i  + M 

/u  s s 
f r o 

a 

)-i 

l1  + 

Jv  rV 

’ a 

-r- 1 

f 

- 

V“v) 

2H  f(s) 
o 


1 + 0 


1 - 


/u  e s\ 


Pc 


2 + od  V— 
*eo 


if  s + 0 , }(  -T—°-  , yj'ucs  d <<  1 


2H  f (s ) 
o 


i - fis. 

" a 


s . , -Zvvas  d 
1 + e 


1 - e 


-2/pos  d 


VU  £ s 

— ^ — <<  1,  \Jvas  d {<  1 


N>|Q 


From  these  expressions,  various  approximate  simple  formulas  can  be  derived 
for  fields  in  the  time  domain,  as  will  be  shown  in  the  following.  Such 
formulas  serve  to  predict  the  approximate  physical  behaviors  of  the  fields 
ar^  to  give  insights  and  provide  cross  checks  to  the  numerical  results. 
Only  response  to  a delta-incidence,  H^110  * 6(t  - J),  is  examined  in  detail 
analytically.  Responses  to  other  incidences  can  be  obtained  by  a convo- 
lution. In  particular,  the  response  to  a "narrow"  incident  pulse  is 
obtained  by  simply  multiplying  the  time-integrated  area  undbr  that  pulse 
by  the  5-response.  The  condition  for  such  a "narrowness"  is 

J.  lies  in  {t^  } and  (t  >>  ltQ)  (12) 

where  {tK  J } ara  the  times  in  which  the  S-response  expression  is  valid 

and  At  is  the  time-width  of  the  incident  pulse, 
o 

2.1.2  Fields  Near  the  Incident  Surface  in  a Thick  Slab 


First,  for  incident  waves  with  frequency  contents  not  too  large  nor  too 
small 


VT 


« V-2— 

Keo 


^ tTTct 


/via"  d 


(13) 


or,  equivalently,  at  times  not  too  early  nor  too  late 


VT „ 


V7#>  V^d 


(13*) 


Equation  (12)  clearly  gives  H(s,0)  * 2Hq  f(s),  a total  magnetic  field  at 
the  incident  side  of  a highly  conductive  thick  slab  being  twice  the 
incident  value,  as  it  should  be.  If  the  incident  pulse  has  its  peak 
time  within  (12*),  as  is  usually  the  case  for  EMP,  then  at  position 
2-0  the  magnetic  field  has  its  peak  time  the  same  as  that  of  the  inci- 
dent one,  but  has  its  peak  value  twice  that  of  the  incident  one. 


Second,  under  a restriction  on  the  frequency  or  time  ranges  that  differs 
slightly  from  (13) , 


/s' 


s » 


d 


or  equivalently 


ur  £o 
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(14) 


(14*) 


and  at  positions  in  the  slab  near  but  not  on  the  left  surface  and  not 
close  to  the  right  surface  such  that 


« z <<  /6  d 


(15) 


the  field  can  be  obtained  from  (9)  by  ignoring  terms  containing 

exp (-/yas  d).  Such  a magnetic  field,  for  a 6-incidence,  is  (Appendix  C) 


H(<5)(t,Z) 


^ 0 


-U02  /4t 


/7Tt 
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(16) 


As  a function  of  time,  at  a given  z satisfying  (15),  the  magnetic  field 
(16)  has  its  peak  value 
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i j at  the  peak-time  (with  given  z as  a parameter) 
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| H Obviously,  the  condition  (15)  on  z ensures  that  satisfies  the 

1 condition  (14*)  on  t,  and  therefore  the  results  (17)  are  valid  under 

1 ||  the  sole  restriction  (15). 

I ~t  Viewed  differently  as  a function  of  position,  at  a given  t satisfying 

i 

— ^ 

■ % 

1 **  (14' ),  the  magnetic  field  (16)  has  its  maximum  value 
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; which  diffuses  to  the  right  with  velocity 
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v(6)  (t)  = ~ z(6),  - cX-  (18c) 

- • z , dt  max/t  *2yot 
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■ 
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Again,  the  condition  (14* ) on  t ensures  that  z j^axyt  satisfies  the 
condition  (15)  on  z,  and  therefore  the  results  (18)  are  valid  under 

the  sole  restriction  (14!). 

LI 

LI 


2.1.3  Fields  at  the  Shielded  Surface  of  the  Slab 


For  a highly  conductive  slab  satisfying  the  first  inequalities  in  (14) 
and  (14’),  without  any  restriction  on  the  slab  thickness,  the  5-response 
transmitted  magnetic  field  at  2 * d is  (Appendix  C) 
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At  times  that  (19)  is  valid  and  converges  fast, 
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(19)  reduces  to 
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This  transmitted  field  has  its  peak  amplitude 
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(1  ~ 28) 
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at  a peak  time 
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yad2S  = pod2 


i - m 


- 2.31  x 10 
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(22b) 


The  time-width  of  the  peak  at  10  strength  is 


(Atpk/d)i/io  ^0.3  yod2  - 7.54  x 10~6 


“'Llo7)^3) 


(22c) 


Notice  that  the  second  inequality  of  (20)  is  always  satisfied  by 
an<*  C^e  fifst  one  is  also  satisfied  by  t^/j  if 


1 « 


(23) 


For  real  shielding  problems,  (23)  is  always  amply  satisfied.  Thus, 
results  (21)  and  (22)  are  valid  approximations  in  the  realistic  time 
interval  of  interest  (20). 

Before  going  into  the  nonlinear  case,  we  make  the  side  remark  that  the 
transmitted  field  H(s,d)  for  a slab  of  thickness  d is  much  smaller  than 
the  transmitted  field  H (s,d)  at  a depth  d in  a semi-infinite  (half-space) 
slab  of  the  same  material.  In  fact,  for  cr  >>  i^c^s,  their  ratio  is 
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Intuitively,  this  is  clearly  plausible  because  the  diffused  field 
reaching  z » d will  leave  the  slab  surface  and  propagate  into  z > d much 
faster  for  the  d~thick  slab  case  than  it  does  for  the  semi-infinite  slab 
case,  and  therefore  the  H(s,d)  has  lass  opportunity  to  pile  up  than  the 
H^(s,d)  does. 

2.2  THE  NOKiXNEAB  u CASS 

For  a ferromagnetic  material,  the  permeability  depends  on  the  magnetic 
field  strength, and  it  is  the  differential  permeability 

g-  0(H)  i »S<H>  „„  (25) 

that  enters  the  field  equation  [5] . In  a strict  sense,  hysteresis  makes 
dB/dH  not  a single  valued  function  of  H.  However,  for  the  transient 
field  behaviors,  not  the  steady  state  behavior,  that  we  are  investigating, 
we  can  use  approximately  [6]  the  magnetization  curve  B(H)  to  get  yR(H). 
Such  a u_ (1!)  for  a typical  iron  saturates  at  a magnetic  field  strength 
„ o£  o*et  rf  ~ _ .<  WM.  . 104  » 

1G3  to  a ur(H>Hc)  V 10  to  1,  in  a range  of  change  in  H of  the  order  of 
tens  of  amp/meter. 

A very  simple  expression  to  approximately  fit  such  a magnetization  curve 
can  be 


VH) 


1 + 


uro  : 1 

1 4-  ea(H'H0 


(26) 


where  yRo  >>  1 and  eaHc  1>  1.  This  fit  gives  a ^(0)  ^ yRQ,  a 
u_(H>>H  ) 'v  l,  and  makes  the  saturation  transition  occur  in  a range 

R C 

AH  v 1/a  about  H H . The  no-saturation  case  is  simply  represented 

c 

by  H -►  «. 
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In  the  following,  we  first  snake  several  theoretical  remarks,  then  solve 
the  nonlinear  y problem  numerically,  and  finally  establish  the  agreement 
between  the  numerical  and  the  analytical  results. 

2.2.1  Theoretical  Remarks 

For  the  linear  y case,  y being  a constant  throughout  the  whole  slab  and 

independent  of  time,  at  fixed  positions  in  the  slab  a smaller  y results 

in  a stronger  diffused  field  («  if’\  see  (17a)  and  (22a))  being  diffused 

to  there  in  a shorter  time  (*  y,  see  (17b)  and  (22b)).  Viewed  differently  at 

fixed  times,  the  spatial  profile  of  the  field  diffuses  and  reaches  its 

"equilibrium”  shape,  peaked  and  symmetric  about  the  center  of  the  slab 

-1/2 

after  the  incident  pulse,  faster  (diffusion  distance  and  velocity  * u , 

see  (18)).  As  a result,  a constant  smaller  y not  only  enhances  the 

diffused  and  transmitted  c. aids,  but  also  enhances  the  higher-frequency 

part  more  than  it  does  the  lower  frequency  part.  The  latter  statement 

1/2 

can  be  seen  for  the  case  of  interest  (y  e s/cr)  <<  1 from  (Appendix  D) 


I H(s>,d) 

\ Sfa'JTd) 

where  y>  > \i<  and  s>  > s<,  or  from  (16),  (17),  (21),  (22)  directly. 

Now  for  a slab  with  a nonlinear  y,  which  saturates  to  smaller  values 
where  the  field  is  stronger,  the  diffusion  results  are  very  different 
but  can  still  be  carefully  extrapolated  from  the  linear  results.  First, 
at  constant  times  the  spatial  profile  has  its  strong-field  center  part 
diffuse  faster  chan  the  low-field  edge  part.  Thus  the  strong-field 
spreads  out  in  a wider  range  and  retains  a lower  value  than  it  does 
without  saturation,  and  it  overtakes  but  is  "confined"  by  the  low-field 
edge  part,  in  a manner  somewhat  similar  to  a shock  phenomenon.  in  short, 
wLg  ssturs t ion  slakes  strong  fisids  diffuse*  s^isily*  -3nd  thus  acts  to 

distribute  the  field  more  evenly  and  mitigates  the  build-up  of  a localized 
strong  field. 


Sh(d/oy<s<)  Shid/cy>s>) 
SfTCdVoy^sp  Sh(d/oy;sp 
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(27) 
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Second,  if  the  pulse  tiae-vidth  is  long  compared  with  the  saturated  peak.-* 

diffusion-time  through  the  slab,  t , with  u *•  u , (see  (22b) 1 

pk/d  * ^saturac ed'  ^ J * 

then  the  diffused  and  transmitted  fields  behave  the  same  as  if  tbs  whole 
slab  has  the  smaller  permeability  ^saturated — wit*1  l&zgst  peaks  and  shorter 
times  as  described  above  in  the  beginning  paragraph.  However,  if  the  pulse 
width  is  short,  i,e,. 


dt  < t . 
o pk/d,u 


i - m 


saturated 


u . ad" 

saturated 


(28) 


and  the  slab  is  thick  relative  to  the  incident  wave  such  that  the 
transmitted  magnetic  field  is  much  smaller  than  the  Hc>  i.e., 
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(29) 


the  maximum  field  at  given  times  and  the  peak  transmitted  field  at  z * d 
becomes  smaller  than  they  are  without  saturation,  due  to  the  first  effect 
just  mentioned.  For  a highly  conductive  slab  that  has  o » us  c y0(B), 
saturated  or  not,  the  amount  of  field  admitted  into  the  slab,  from  (11"), 
is  virtually  independent  of  y.  The  maximum  field  at  a given  time  is  roughly 
inversely  proportional  to  its  spatial  spread  at  that  time.  Thus,  approxi- 
mately we  have  the  ratio  of  the  maximum  transmitted  fields  with  saturation 
to  that  without  saturation  (Appendix  E) 
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(30) 


where  the  superscripts  N.L.  sat.  and  no  sat.  denote,  respectively,  the 

case  of  nonlinear  saturation  and  no  saturation,  the  t is  the  time  at 

c 


which  the  staxisuas  magnetic  field  for  the  no-saturation  case  decreases 
through  the  8^,  and  zmr.vr- 1 is  the  distance  the  maximum  field  reaches  at 
time  t . From  (2-1$)  and  of  course  under  Its  validity  condition  (14'), 
re  <-an  express 


(no  sac.)  „ 
2max/tc 


llilRl  • &t  \l/2 


ltcvo 


o'  8, 


(31) 


The  value  of  \ larger  than  z^/^t  * as  discussed  in  the  first 

remark,  cannot  be  obtained  analytically  in  the  present  analysis,  but 
its  numerical  value  can  be  used,  together  with  (31),  in  (30)  to  relate 
the  nonlinesar  maximum  transmitted  field  to  the  linear  one. 


Third,  the  and  the  (dt^ k/d>l/l0  ^or  C^e  nonllJie3r  case  is  about  the 

same  as  that  of  the  no-saturation  case,  as  long  as  (23)  and  (29)  are 
satisfiea.  For  under  such  conditions  the  field  is  well  dispersed  below 
H long  before  it  diffuses  to  z » d,  and  thus  it  is  the  uasaturated  y_ 

w SO 

that  controls  the  tpk/d* 


Finally,  for  the  EMP  shielding  cases  of  practical  interest,  (28)  and 
(29)  are  satisfied.  This  is  easily  seen  by  substituting  typical  numbers 
into  chose  expressions.  Thus,  the  above  observations  are  practically 
applicable. 


2.2.2  Numerical  Method  and  Results 


Tba  numerical  code  DIFI2SN  solves  the  one-dimensional  nonlinear  y diffusion 
problem  by  finitely  differencing  (5)  and  (?)  with  yR(H)  replacing  the 
constant  yR.  An  implicit  "T"  finite  difference  scheme,  stable  in  the 
round-off  error  and  the  differencing  grid  sizes,  is  used  [?].  The  code  takes 
as  inputs  any  nonlinear  function  uD (H) , any  incident  pulse  shape  H f(t), 
and  the  properties  zQ  or  the  ambient  medium  and  c,  o,  d of  the  slab. 
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As  outputs  it  gives  at  each  glsae  step  the  fields  at  some  selected  fixed 
positions,  including  of  course  2*0  and  z a d,  the  value  and  the  location 
of  the  maximum  j&gguetic  field  (in  0 < 2 < d)  at  that  lisas,  and  tits  two 
locations  between  vhieh  the  magnetic  field  at  that  tine  exceeds  a selected 
value  such  as  the  Hc  about  which  the  saturation  occurs,  A listing  of 
the  code  01FUSN  is  included  in  Appendix  F. 

numerical  results  for  a number  of  parameter  values  of  practical  interest 
ere  obtained  and  plotted.  In  these  results,  the  p^(H)  of  (26)  is  taken 
to  represent  the  non-linear  permeability,  and  the  expression 

HiflC(t,z)  - H sin2  u (t— ).  0 < t-jr  < — 

G O C — C — aJ 

» 0,  otherwise  (32) 

is  used  as  the  incident  wave.  The  resulting  plots  and  their  comparisons 
with  analysis  are  given  in  the  following. 

The  results  for  a typical  example  of  iron  shielding  with  a pR(H)  whose 

y_  « 10\  d * 400  asap/m,  and  a * 1/50  m/amp,  a thickness  d » 3 mm, 
xo  c , - 

a conductivity  <3  •*  10  mho/m,  and  an  incident  H “10  aap/m  and 

t 0 

» 3 x 10  rad/sec  are  plotted  in  Figure  2 to  Figure  5.  Corresponding 

results  for  this  same  case  but  with  constant  y„  » * uR(0)  ar.d 

uf  » 1 * ur(*)»  limiting  values  of  the  nonlinear  u^(H),  are  also 

computed  and  plotted  together.  The  plots  show  the  maximum  magnetic  field 

H(t,z  , ) of  the  diffusion  profile  first  increases  as  twice  the  incident 
MX/  c 

wave  (Figure  2-1)  at  the  incidenc  surface  of  the  slab  (Figure  2-2),  then 
breaks  away  from  the  incident  wave  and  decreases  (Figures  2-1,  3-1,  4-1) 
while  diffusing  into  and  toward  the  center  of  the  slab  (Figures  2-2,  3-2, 
4-2).  For  the  construe  u cases,  the  values  and  the  location  of  H(t,z  . ) 

nisix/  t 

agree  very  well  with  those  given  in  (IS)  (in  its  region  of  validity  (14  * ) 
of  course)  from  the  previous  analysis  (labeled  curved  in  Figures  2-1,  3-1, 
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H (t,  *MAX  t*  <8fnp/fT>eter) 


itr6  to'5  io-4  io-3 

t (sec) 

Figure  3-1.  The  Maximum  Fields  in  the  Slab  H(t,  z„  .*)  as  Functions 

rnax  I t 

of  Time  (Continuation  of  Fig.  2-1) 
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The  value  of  for  th*  nonlinear  y case  must  be  less  than  the 

larger  of  those  for  the  corresponding  limiting  constant  y cases;  it 
approaches  that  for  **  1 if  the  pulse  has  long  duration  and  strong 
amplitude  to  saturate  the  whole  thin  slab  (i.e.,  (28)  and  (29)  inequalities 
reversed),  and  vice  versa.  But  it  does  not  do  so  mono conically,  because 
of  the  competing  saturation  effect,  causing  the  strong  field  near  the 
z - 0 end  to  diffuse  faster  both  in  the  +z  direction  penetrating  the  slab 
and  in  the  ~z  direction  escaping  the  slab.  The  location  zmavfr  of 
H(ttzma T/i.')  f°r  ^nlinear  U case,  however,  always  lies  between  those 

for  the  corresponding  constant  y cases  (Figures  2-2,  3-2,  4-2),  as  it 
should,  because  the  smaller  saturated  y under  the  peak  permits  it  to 
diffuse  faster.  Further,  the  spatial  region  within  which  the  H(t,z) 
exceeds  the  Hc  of  the  nonlinear  y case  is  plotted  (Figures  2-3,  3-3). 

This  region  is  roughly  the  extent  within  which  the  nonlinear  saturation 
occurs,  and  its  disappearance  marks  approximately  the  end  of  the  satura- 
tion effect. 


The  transmitted  magnetic  fields  H(t,d)  for  various  parameters  are  shown 
in  Figure  5 to  Figure  9.  First,  the  limiting  constant  y results  exhibit 
excellent  agreement  with  the  analytical  formulas  (21)  and  (22) , for  the 
typical  example  (Figure  5)  and  for  other  variations  of  parameters  (not 
plotted),  when  the  (practical)  thick  plate  condition  (23),  thus  the  condition 
(20),  and  the  short  pulse  condition  (28)  are  satisfied.  Second,  satisfying 
the  additional  but  still  practical  condition  (29),  the  H(c,d)  in  the 
typical  nonlinear  example  has  the  same  time  shape  as  but  is  slightly 
lowered  near  its  peak  by  a factor  ''-0.8  from  the  H(t,d)  in  the  same  typical 
example  but  with  a constant  yr  * yRo  (Figure  5).  This  agrees  with  (30) 
as  it  should.  Third,  under  the  restrictive  but  practical  conditions 
(23),  (28)  and  (29),  the  H(t,d)/HQ  with  different  Hq  substituted  in  the 
typical  nonlinear  problem  (Figure  6)  decreases  slightly  near  its  peak  for 
larger  Hq,  approximately  according  to  (30),  but  has  its  time  shape 
virtually  unchanged.  This  is  as  expected  from  Section  2.2.1.  Similar 
variations  in  w , in  view  of  the  multiplicative  factor  vHQ/ (2wq)  to 
convert  the  5-response  to  the  narrow-pulse-response  (32),  expectedly 
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give  similar  results  (Figure  ?).  Fourth*  under  these  ease  conditions  but 

varying  the  slab  thickness  d (Figure  8)  and  the  = y^(0)  (Figure  9), 

respectively,  in  the  typical  nonlinear  example  shows  a time  scaling 

2 

proportional  to  d , and  an  amplitude  dependence  approximately  as 


"oy 


~1„3 
Ro" 


but  slightly  less  than  such  a scaled  value.  Again,  this  is 


expected  from  (22)  and  Section  2.2.1. 


Me  conclude  this  numerical  section  by  stating  again  that  for  short  pulses 
and  thick  slabs,  under  condition  (23),  (28),  and  (29),  the  transmitted 
magnetic  field  for  the  nonlinear  u obeys  approximately  (22) , with  * 
u^,  but  becomes  slightly  smaller  as  corrected  by  (30). 


SECTION  3.  THE  CYLINDRICAL- INCIDENCE  SLAB  PROBLEM 


Consider  Figure  10,  an  incident  cylindrical  TEM  wave 

- s?  £(‘  - 1)  «"■> 

E>-o-s) (33b) 

o 

carrying  a z-flowing  current  IQf(t  - z/c)  on  the  surface  of  a perfectly 
conducting  wire  of  radius  a impinging  upon  a slab  as  shown.  In  the 
cylindrical  coordinates  (p,<j>,z),  the  <£— symmetry  makes  fields  functions 
of  (t.p.z).  The  parameters  in  Figure  10  have  the  same  meanings  as  in 
Figure  1. 

Before  solving  the  problem,  we  must  make  several  remarks  here.  First, 
one  can,  of  course,  solve  the  problem  numerically  either  by  finitely 
differencing  the  field  equations  subject  to  the  bounds  onditions  or 
by  finitely  patching  the  perfect-conductor  mixed  surface  integral  equations 
to  solve  for  the  surface  current  density  on  the  wire  and  then  to  obtain 
the  field  from  the  surface  current  density.  Second,  with  its  end  surface 
present  the  semi-infinite  cylinder  does  not  have  its  surface  as  one  of 
the  coordinate  surfaces  of  the  eleven  coordinate  systems  that  permit  the 
separation  of  variables  for  the  Helmholtz  equation.  Thus  in  using  the 
familiar  method  of  solving  by  the  separation  of  variables  and  summation 
of  the  products  of  the  eigen  functions  (in  thi3  case  the  Hankel  transform 
for  o or  the  (2-sided)  Laplace  transform  for  z) , not  only  are  the  coeffi- 
cients mixed  by  the  boundary  conditions  but  also  the  expansion  does  not 
converge  at  the  edge  of  the  cylindrical  end  where  the  electric  field  goes 
to  infinity  (although  weakly  if  edge  condition  of  finite  energy  is  imposed) . 
But  to  match  boundary  conditions  means  that  we  would  have  to  manipulate 
and  evaluate  the  coefficients  of  a series  at  precisely  the  place  where  the 
series  representation  is  not  valid.  As  a result,  no  con.  istent  solution 
can  possibly  be  obtained  by  such  a familiar  method.  The  rigorous  analytical 
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solution,  even  to  the  simplified  problem  without  the  presence  of  the  slab, 
is  therefore  still  an  unanswered  question  [8]. 

In  the  following,  we  solve  the  cylindrical-incidence  problem  by  an 
approximate  method  whose  validity  range  is  established  to  be  of  practical 
interest.  The  resulting  formulas,  being  simply  related  to  those  for  the 
one-dimension  problem,  enable  us  to  make  use  of  results  in  Section  2 for 
the  present  cylindrical  problem. 

3.1  ANALYSIS  FOR  THE  CONSTANT  y CASE 


To  find  the  diffused  fields  and  the  transmitted  fields  in  z > 0,  the  only 
physical  quantity  we  really  need  from  the  incident  z < 0 side  is  the 
tangential  field,  either  electric  or  magnetic,  at  z * 0.  Although  it 
does  not  matter  which  one  we  start  with  if  we  can  solve  the  problem  exactly, 
it  does  make  a difference  in  solving  the  problem  approximately  depending 
on  which  approximation  can  be  made  more  easily  and  more  accurately. 

Now,  for  the  real  problem  both  the  "wire  region"  p < a;  z < d and  the 
"surrounding  slab  region"  p _>  a,  0<z<d  are  highly  conducting.  The 
magnetic  field  at  the  "inside"  boundary  p >'a-6  , z = 0,  where  6^  is  the 
skin  depth  of  the  wire  region  and  is  <<  a,  is  clearly  much  smaller  than 
the  at  the  "outside"  boundary  p > a,  z ■ 0,  regardless  of  the  existence 
of  the  conducting  surrounding  slab.  In  fact,  without  the  slab's  presence 
this  has  been  used  to  justify  the  incident  current  I being  associated 
to  a cylindrical  TEM  wave  outside  the  wire,  and  with  the  slab's  presence, 
and  reflection,  this  is  even  more  pronounced.  Moreover,  the  area  of  the 
"inside"  boundary,  , is  much  smaller  than  the  area  of  the  "outside" 
boundary,  ».  This  makes  the  diffused  and  transmitted  fields,  being  the 
surface  integration  of  the  cross  product  of  the  He  and  the  gradient  of 
the  Green's  function,  depend  even  less  on  the  "inside"  magnetic  field. 

Thus,  the  H.  at  the  "outside"  boundary  alone  predominatly  determines  the 

v 

diffused  and  transmitted  fields.  However,  a similar  result  does  not  hold 
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for  the  boundary  electric  field-  Because  the  "inside"  boundary  electric 
field  E , though  negligibly  small  relative  to  the  "outside"  boundary  one 
when  there  exists  no  conducting  surrounding  slab,  is  not  comparatively 
negligibly  small  when  the  conducting  surrounding  slab  is  present  and  makes 
the  "outside"  boundary  Ep  also  vanishingly  small.  Consequently,  we 
must  use  the  magnetic  matching,  not  the  electric  matching,  if  we  use  the 
"outside"  boundary  fields  only  as  an  approximation. 

Finally,  for  the  case  of  narrow  incident  pulse  and  thick  slab  of  practical 
interest  ((28)  & (29)),  the  approximate  "outside"  boundary  magnetic  field 

can  be  found  easily.  Since  the  slab  makes  the  H,  at  p > a,  z * 0 not 

$ 

sensitive  to  the  geometry  beyond  z > d,  the  there  is  virtually  unchanged 
if  we  extend  the  wire  beyond  z ■ d to  z » «.  Furthermore,  since  the  wire 
is  itself  highly  conducting,  as  far  as  fields  on  it  are  concerned  we  can 
replace  it  by  a perfectly  conducting  wire.  The  following  analysis  follows 
such  an  approximation  procedure. 

3.1.1  Formulation  and  Analysis 

If  the  perfectly  conducting  wire  extends  to  all  z,  then  the  incident  TEM 

wave  (33)  produces  only  TEM  waves.  Their  wave  forms  for  a constant  p slab 

are  just  the  one-dimensional  results  with  H reolaced  by  H a/p  where 

o • o 

Hq  - Io/(27ra),  and  with  Ex  and  substituted  by  E q and  H^,  respectively. 

In  particular,  from  (11) 


^(<*>  wire) 


(s.p ,0) 


JiL  . 2f (s)[(l+a)  - (l-a)e’2/^  dj  > a 
2lrp  (1+a)2  - (l-a)2e~2,/ir5i  d 


= 0 , p < a 


(34) 


where  a = (y  e s/a) 
r o 


1/2 
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Sow,  for  the  real  problem  of  a terminated  wire  with  ^-symmetry , the  TM 
wave  cannot  exist  in  the  region  z > 0 because  of  its  necessarily  accompanying 
singularity  at  P * 0.  Further,  a TE  wave,  under  the  condition  3/3$  = 0, 
has  only  , and  H2  as  the  non-vanishing  components.  But  the 

boundary  conditions,  the  continuity  of  tangential  fields  at  z * 0 and  z « d, 
do  not  couple  this  TE  mode  to  the  incident  TEM  wave  at  all.  Thus  the  TE 
wave  has  its  existence  independent  of  the  incident  wave  and  can  be  taken 
as  identically  zero.  Finally,  a TM  wave  does  exist,  with  non-vanishing 
components  H?1*,  E?*,  and  eP*‘  In  z > 0 , this  is  the  only  mode  existing,  and 

T P Z 

its  superscrip  s are  omitted  in  the  following. 


Similar  to  the  one-dimension  problem,  the  diffused  fields  in  the  highly 
conducting  slab  0 < z < d satisfy  (2)  and  obey  the  equations 
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(35b) 


(35c) 


In  z > d they  obey  similar  equations  as  (35)  but  with  a replaced  by 
eQ3/3C.  To  solve  the  problem,  in  addition  to  the  Laplace  transform  in  time 
we  make  use  of  the  Hankel  transform  in  p [9]. 

<Kp)  “ / dK  J (Kp)<j;(K),  p > 0 (36) 

Jq  v 

Now,  requiring  as  usual  a + z traveling  wave  in  z > d,  finite  fields  at 

* * (M 

p * 0,  continuous  E and  H,  at  z = d,  and  that  H,  (s,p,0)  = Hv  *"  /(s,p,0) 
P $ 9 

of  (34),  the  magnetic  field  in  0 < z < d is  (Appendix  G) 
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where 


H (S,P,z) 


<*K  J^Kp) 


?<K,s)e 


-K>z 


+ n(K,s)e 


K>z 


(37) 


/(l+8)eK>d  \ 

q(K,S>  \ j / 

\-(l-g)e  K>d/ 
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(38) 
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IT 
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Real{ K>)  > 0 , 

(39c) 

/ 2 2 
2 /K  +u  e s 

o o * 

Real{  K<)  > 0 

(39d) 

The  electric  fields  in  0 < z < d are  given  immediately  by  (37)  and  (35a,b) 
The  fields  in  z > d are  (Appendix  G) 


H (s,p,z)  - se  / dK  Kp(K,s)JOCp)e”Ki<Z"d) 
9 °Jq  1 

(40a) 

E (s,p,z)  -/  dK  KK<P(K,s)J.  (Kp)e“K:<(2_d) 

P Sq  1 

(40b) 

/*“ 

E (s*p,z)  - / dK  K2p(K,s)J  (Kp)e‘K<<Z'd) 
* JQ  o 

(40c) 

where 


p(K,s) 


q(K>s) 

se  K 
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eK>d(l-bs)  - e"Kid(l-S) 


(41) 
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For  a highly  conducting  (a  <<  1,  condition  (13))  thick  (od/u  '/e  1 

condition  (23)>3lafa,  £,  n,  and  p reduce  to  (Appendix  G) 
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We  examine  next  such  cases  in  detail. 


From  (40b),  (41),  and  (43),  the  transmitted  radial  electric  field  at 
z * d is  (Appendix  H) 


and  becomes 
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for  a 5-incidence  H^Ct.p.z)  » 6<t-z/c)a/P.  Here  9 (vjt)  is  one  of  the 
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Theta  functions  [10] . 


(45a) 


(45b) 


(45c) 


Similarly,  from  (40a)  the  transmitted  magnetic  field  is 
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rr~ 

(<t,p ,d)  Ep(t,p,d) 


+ H (t.P.d) 

9,  dev 


(46) 


Here,  ^^(t.p.d)  is  the  part  of  magnetic  field  which  deviates  (in  its 

relation  to  E ) from  that  of  a pure  TEM  wave  (Appendix  H) 
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(43a) 
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Thus  at  large  p and  small  p the  total  magnetic  fields,  respectively 


are 


2 

UO0  } 


[*  visiVf  4--V7-  . is 


if  p » a,  e(ii0f>‘t)/4  » 1, 


(49a) 
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if  p « a (49b) 

for  a Hinc(t,p,z)  - <5(t  - |)  a/p. 

Comparing  (45)  to  (49)  with  the  one-dimension  results  (4)  and  (19),  we 
see  that 


H<5)<t,p,d)  ■v  (t.d)  * [l  - e(-yop2)/  (4t) 

$ one-D  p l 


if  p » a,  e(uoc  C)M  » 1 
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and 


E(<S)(t,p,d)  • E(5)  (t,d)  - 

one-D  p 
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The  relations  (50a)  and  (51a)  for  large  p a are  intuitively  plausible: 
at  early  times  before  the  effect  of  the  wire  region  diffuses  there, 
t <<  yap  /4,  the  fields  do  not  feel  the  termination  of  the  wire  and  are 
just  the  diffused  cylindrical  TEM  wave  with  the  one-dimensional  fields 
multiplied  by  a/p,  and  at  later  time,  the  fields  become  smaller  than  such 
a TEM  version  because  they  diffuse  into  the  wire  region.  The  normalized 
ratios  H^<5)(t,p,d)/(H^_D(c,d)  a/p)  and  Ep(5)  <t,p  }d)/(E^_D(c,d)  a/p) 
as  functions  of  a normalized  time  X =*  c/(yop^)  show  (Figure  11)  such 
behaviors.  The  one-dimensional  peak  diffused-through  time  “ 

t^/d/(yop2)  (see  (22b))  as  a function  of  d/p  shows  (Figure  11)  that  for 
d/p  ^ 1 the  diffused  fields  are  virtually  a pure  cylindrical  TEM  wave  up 
to  the  peak. 


For  small  p <<  a.  the  H^(t,p,d)  and  E^(t,p,d)  given  respectively  fay 

9 P 

(50b)  and  (51b)  go  to  zero  at  least  as  fast  as  p , a consistent  behavior 

with  a finite  Ez  and  a finite  total  (conducting  and  displacement)  current 

density  near  p ■ 0,  and,  of  course,  start  from  zero  at  t * 0.  The  ratios 

(t,p  »d)/Hpp^_D(t,d)  and  (t  ,p  ,d) /Epp^_D(t  ,d)  show  such  behaviors 

(Figures  12-1,  12-2,  plotted  similarly  to  Figure  11  but  with  p replaced 

2 

by  a for  the  normalized  time  t/(yc?a  ) and  with  p as  a parameter). 
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RATIOS  OF  FIELDS 


1 


X a t /(pap2) 


Figure  11.  The  Nonullzed  Ratios  p,  d)/(H®  (t,  d)  |)  and 

EP  (t,  p,  d)/(E1>i0{t,  d)|)from  thep»a  Expressions  (50-a) 
and  ( 5 1 - a ) as  Functions  of  Normalized  Time  t/{pap2),  and 
the  Normalized  1-0  Peak  Time  Xpk(<j  * t ^/(pop2)’ 
as  a Function  of  (d/p) 


/a  * 3 x 10' 


P/a  <•  10* 


p/a  - 3 x 10'' 


P/a  - 10 


pk|d  - 9.13  x 10'2  {d/ai2  from  (22b) 


Y = t/ipaa*) 


Figure  12-1.  The  Ratios  H^(t,  p,  d)/Hj_0(t,  d)  from  the  p«a 
Expression  (50b)  as  Functions  of  Normalized  Time  t/(paa^) 

and  the  Normalized  1-0  Peak  Time  Ypkjd  = tp(cj d/(pca2) 
as  a Function  of  (d/a) 


\ P/a  * 3 x 10' 


J'  \ \?; 0 \ 


.Ypk|d  = 9.18  x 10'2(d/a)2 
FROM  (22b) 


Y = tApcra*) 


Figure  12-2,  The  Ratio  Ep'(t,  p,  d)/£^Q(t,  d)  from  thep«a 

Expression  (51b)  as  Functions  of  the  Normalized  Time  t / 

2 

{•p'ja  ),  and  the  Normalized  1-D  Peak  Time  Y^^  j ^ = 
tpk|  as  a Function  of  (d/a) 


Now  consider  p v a.  From  (45b),  the  £v  'l(t,p>>a,d)  expression  is  also 

valid  for  E^  ;(t<<yaa  ,pva,d).  Similarly  from  (45c),  the  expression 

(t,p<<a,d)  is  valid  for  (t>>yoa^,p'\«a,d) . Thus,  for  p **  a,  the 

ratio  ' (t,a,d)/Egn^_jj(t,d)  at  early  times  t <<  \ioa  is  given  by  the 

large  p curve  (the  F-curve,  X « 1 part  in  Figure  11  with  p * a)  and  at 
2 

late  times  t >>  yea  is  given  by  the  small  p curve  (the  Y » 1 part  of 

2 

the  p * a curve  in  Figure  12-2).  At  t v y a a , they  do  join  smoothly 
at  M).2  (the  E-curve  at  X 1 in  Figure  11  and  the  p/a  * 1 curve  at 
Y v 1 in  Figure  12-2)  which  should  be  the  right  value.  Similar  conclu- 
sions hold  for  the  (t,pva,d)/H<‘<5^  n(t,d).  This  is  because  at  P ^ a 

<?  one-u 

the  part  of  magnetic  field  which  diverts  (see  (46))  from  the  TFM-like 
behavior  is  negligible  for  early  times  and  for  late  times,  respectively, 
in  the  large  p expression  (50a)  and  in  the  small  p expression  (50b). 


We  conclude  by  restating  that  (50)  and  (51)  (plotted  in  Figures  11,12) 
are  valid  relations  between  the  cylindrical  and  the  one-dimensional 
solutions  for  highly  conducting  thick  slabs  of  practical  interest.  From 
these  relations,  we  can  multiply  the  one-dimension  fields  obtained  in 
Section  2 by  a cylindrical-effect  factor  to  make  predictions  for  the 
cylindrical-incidence  problem. 

3.2  THE  NONLINEAR  y CASE 

The  nonlinear  saturation  effect  can  only  occur  in  a region 


a v p £ — a 


From  Section  2.2  we  have  found  that  for  the  practically  interesting 
cases  of  highly  conducting  thick  slabs  with  narrow  incident  pulses 
(conditions  (23),  (28),  (29)),  saturation  only  reduces  slightly  the 
transmitted  fields  but  leaves  intact  their  time  shapes.  Such  are  the 
affects  the  cylindrical  transmitted  fields  in  region  (52)  will  experience. 
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with  (50)  and  (51)  still  being  the  formulas  to  relate  to  those  of  the 
one-dimension  nonlinear  problem  which  are  solved  semi-analytically  or 
numerically  - 

For  p not  in  the  region  (52),  the  constant-p  results  in  Section  3.1  apply 
directiy. 


Finally,  we  shall  conclude  with  an  example.  Consider  the  typical  problem 

in  Section  2.2.2,  but  with  a wire  of  radius  ,vQ.5  cm  carrying  an  incident 

3 

1Q  * aHQ  ■ 3-14  x 10  amp.  The  transmitted  magnetic  field  in  the 
region  p ^ 1.25  m is  approximately  the  «*  103  amp/m  curve  in  Figure  6 
multiplied  by  the  cylindrical  factors  (50)  or  its  plots  (Figures  11,  12-1). 
For  example,  at  p a,  the  transmitted  field  has  a peak  value 


Vtpk/d’p'V'a’d)  ^ 10' 


x 3.8  x 
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11 


with  a peak  diffusion  time  'vl.l  x 10  ^ sec  and  a diffused  pulse  width 
(ht)1/iQ  4.0  x 10  1 sec,  same  as  for  the  one-dimensional  case.  For 
p > a,  the  peak  field  is  smaller  by  a factor  ''Wo,  but  that  factor  approaches 
MO. 8)  ^a/p  as  p >>  1.25  m.  At  p < a,  the  peak  decreases  from  (53)  at 
P 'v  a to  zero  at  p *>  o,  roughly  proportional  to  p as  p £ a and  according 
to  (50b)  as  p<<a. 
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APPENDIX  A.  BOUNDARY  CONDITIONS  (7) 

From  the  Maxwell  equations,  by  ignoring  e 3/3t  compared  with  a,  we  have 
in  the  slab  0 < z < d 


3E 


3H 


3z  ~vr  po  3t 


(A-l) 


crE 


3K 

3z 


(A~2) 


which  immediately  give  (5).  Now  the  boundary  conditions  at  z - 0 


are 


H(t,o)  - H [f (t)  + g(t)] 


(A- 3) 


-V5 


3^  H(t,o)  - H0o[f(t)  - g(t)] 


(A-4) 


which  immediately  give  (7a) , and  the  boundary  conditions  at  z * d are 


H(t,d)  * HoT(t) 


(A-5) 


3z 


H(t,d) 


-# 


Ho0T(t) 


(A-6) 


which  immediately  give  (7b) 


0 
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APPENDIX  B.  DERIVATION  OF  FIELDS  IN  THE  SLAB,  ONE-DIMENSIONAL  PROBLEM 


In  che  Laplace  domain,  (5)  reads 


— r H(s,z)  - nos  H(s,z)  ” 0 


(B-l) 


I . 


where  Q(c<0,z)  = 0 has  been  used,  and  (7)  reads 


fj-  H(s»o)  - oH(s, 


o)  - -2y^Ho<rf(s) 


(B-2) 


3z 


H(s,d)  + 


d)  * 0 


(B-3) 


Thus,  the  solution  of  (B-l)  is 


H(s 


, z)  - A(s)  z + B(s)  e1^  s (B-4) 


where  A(s)  and  B(s)  are  to  be  determined.  Now,  from  (B-2)  we  have 


(-A  + B) 


yj  uos  - y-^- 


H f | 
f e o 


a (A  + B)  - -2cr  f (s)  (B-5) 

o o 


and  from  (B-3)  we  have 


<-a  * ♦ b ♦ V5  0(4  4 


+ B e/u0S  d)  - 0 


(B-6) 


Solving  for  A(s)  and  B(s)  from  (B-5)  and  (B-6),  substituting  their  values 
into  (B-4),  and  using  (6)  yields  (9)  and  (10). 
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APPENDIX  C.  DERIVATION  OF  3 CUE  APPROXIMATE  FORMULAS 


: S 

» a 


First,  under  the  condition  exp [ /yes  (z-2d)]  <<  1,  the  magnetic  field 
obtained  by  using  only  the  A(s)  term  in  (9a)  for  a 5-incidence  is 


t i 


S(6)(t,z)  * 2 


V?F 

T r o 


A2'44  A+ti  , 

- a e erfc 


4.1 


(C-l) 


which  reduces  to  (16)  if  2cz  Vy^Te^  >>  1.  Here  k = •'pa  2 and 
a 5 /o/(ur£0)  . 


I 4i 


Second,  from  (11*)  the  result 


H(5)(t,d)  » L*1 


4 


y e s 
r o 


Sh(/ycs  d) j 


;VVlt  [/sh94 (°l‘"2)l 


^r^o  d 

1.  V 

' o dt 

/yo  d * 

I MCfd 


-(n-1/2)  2uod2/t 


i-r  TTjE"2"'1^2-2'1  a 


-yod2(2n-l)2 
4t 


(C-2) 


is  used  in  (19) 
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mstis-i 


it 

Ii 


li 


t i 

i i 

u 


U 

11 


ii 

4 i 


I 5 

4 ♦ 


APPENDIX  D.  HIGH  FREQUENCY  TRANSMITTED  FIELD  RELATIVE  INCREASE 

From  (11*) , for  y * p>  the  ratio  of  the  transmitted  magnetic  field 
at  a frequency  s * s>  to  that  at  s * s<  < s>  is 

^H(s>,d)^  (h^)/(Sh(h^) 

R(y>}  5 \H(s<Jd)/u>  * (M^)/(Sh(h/^)  0 

where  h = d^s“ . Now  we  want  to  prove  that  for  a smaller  y * p<  the 
H(u<)  becomes  larger,  i.e. 


R(u<)  > R(u>) 


(D~2) 


or  rewriting  it  equivalently 


Sh&x  Sh  x 

Sh(£rx)  Sh(rx) 


(D-3) 


Here  x = h/y'Ja”  > 0,  r = /s^/s^  > 1 and  i = ^'y</y>  < 1.  But  the  function 
Sh(2.x)/Sh(irx) , viewed  as  a function  of  £,  is  monotonically  decreasing 
since  r th  y > th(ry)  for  r > 1 and  y > 0.  This  establishes  (D-3)  and 
completes  the  proof. 
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APPENDIX  E-  ARGUMENT  FOR  EQUATION  (30) 


11 

i I 
; * 

i i 

11 


Since  the  slab  is  highly  conducting  with  a » whether  it  be  saturated 

or  not,  the  "amount"  of  field  admitted  into  the  slab  is  about  the  same  for 
either  case  (see  (11*)}.  Now  think  of  the  H field  as  the  density  (one- 
dimensional in  z)  of  some  particles  that  diffuse  in  the  slab.  Then  the 
total  number  of  them 


/ 


Hdz  'v 


<«***> <**> 


* H(t  ,2 


c’  max/ t J max/t 


) z. 


is  about  the  same  without  and  with  saturation,  if  the  pulse  is  narrow  and 
the  slab  is  thick  (see  (28)  and  (29)).  After  both  cases  diffuse  in 
the  same  manner  and  should  roughly  give  a peak  transmitted  field  propor- 
tional to  H(t  , z . ).  This  makes  (30)  plausible. 
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APPENDIX  F.  THE  CODE  DIFUSN 
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APPENDIX  G.  DERIVATION  OF  FIELDS,  CYLINDRICAL  PROBLEM 

The  ^-independent  cylindrical  TM  mode  in  the  highly  conductive  slab 
0 < z < d,  from  (35),  assumes  the  form 


-TM 

a,  » hr1 

<p  <p 


r «Ji(S.)[e(^.)  „<*,.)  e'k2+“M  * 


<G-1) 


E - E™  =~-  |-  H 
p p a 3z  41  a 


- f d K.V?"  + pas  J^(Kp)|^(K, 


s)  e 


->'fc2+yas  z 


+ n(K,s)  e 


*/K2+uos  z 1 


(G-2) 


;TM.  119 


dK  KJ 
u 


0(Kp)  [?  (K,s) 


uas  z 


+ n(K,s)  e 


(G-3) 


where  the  finiteness  of  the  fields  at  p = 0 has  been  used  to  choose  the 
Bessel  function  of  first  kind.  Similarly,  the  fields  in  z > d can  be 
expressed  as 


K = sVo)  ^<s>p.z> 


3 “ 2 

— r - s 

3z- 


r 


dK  • p(K,s)  Jn  (Kp)  e 


-^K^+Po^s2  (z-d) 
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(G-4) 
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K2  p<K,s)  J (Kp)  e *'K2+'*ioc0s2  (z-d) 


% = flJ”  “ ~ S£o  *(s‘p’z> 


3 A, 


se_ 


/*"  dK  • K p<K,s)  JL(KP)  e-/K2+y°C°s2  (z"d> 
•rn 


(G-5) 


E « E^  » _J! v(s  o z) 

P 0 3p3z  t‘S’p,Z} 


J dK  * kVk2  + uoeo  p(K,s)  J} 


(Kp)  e 


Vk2+vi0e0s^  (z-d) 


(G-6) 

where  the  V(a,p,z)  is  the  cylindrical  TM  Hertz  potential,  and  the  require- 
ment that  the  wave  be  positive-z  traveling  has  been  used  to  choose  the 
negative  exponential  expression.  In  the  above,  the  choices  of  branch  cuts  (as 
a function  of  s)  are  such  that 


Real 


Real 


{ Vk“  + uos  } > 0 

(G-7) 

(Vk2  + u c s2  } > 0 
0 o 

(G-8) 
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Now  at  2 * d,  the  continuity  of  H gives 


5(K,a)  e”K>d  + n(K,s)  eK>d  » seQK  p(K,s) 


(G-9) 


and  the  continuity  of  gives 

P 


— K 

~ C-C(K,s)  e“K>d  + n(K,s)  eK>d]  * KK<  p(K,s) 


(G-10) 


where  2 ^C^+uos  and  K,,  = \HL*+uq£o8*  . At  z **  0,  requiring  that  the 
total  magnetic  field  be  the  approximate  value  (34)  and  using  the  integration 
expression 


f iK  h 


(Kp)  Jc(Ka)  - ~ , p > a 


(G— 11) 


* 0 , p < a 


(G-12) 


gives 


C(K,s)  + n(K,s) 
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= q (K,s) 


(G-13) 


Solving  (G-9),  (G-10)  and  (G-13)  gives 


(1  + 8)  e ° 
■a  - g)  e~K>d 


q(K,s) 
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(G-14) 
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p(K,s) 


28 


(G-15) 
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eK>d(l  + 8)  - e^>d(l  - g) 


where  g 5 seoK>/(oK<),  These  give  (38),  (35)  and  (41). 

If  « = (ureQs/a)  ' <<  1,  Chns  B « 1,  and  od  yaQ/e  » 1,  a 

practical  condition  as  stated  in  (23)  which  implies 
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where 


1 £ 


X(K,s)  = 


Vo(Ka)  f(s) 


These  give  (4  2)  and  (43). 


(G-20) 
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APPENDIX  H.  APPROXIMATE  EXPRESSIONS  FOR  CYLINDRICAL  PROBLEM 


it 
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i * 
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h i 


From  (40b)  and  (43) , we  have 
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Thus  the  <S-rasponse  for  an  incident  H * a $(t  - z/c)/p,  i.e. , H0f(t)  * 
5(t),  is 


^s>  <'•».“> 


&•  e401^)]  /"  “ JilKp')  j°(k*> 


(H-4) 


,2a  d . 
<jd  dt  4 


if  p » ain 


(0|^)]  ‘ *■ 
■(••  V£) 


(-yop  )/(4t) 


i . ,2a. 

<L  9 /( 

jtiTtt  \ 

cd 

dt  84\( 

3 ,2/ 
ycd  >1 

if  p « max  (a, 

Similarly,  from  (40a)  and  (43),  we  have 

H^(t,p,d)  - L"1  f dK  Jx(Kp)  J0(! 

^ A 


H^(t,p,d)  » L 


dK  Jx(Kp)  J0(Ka)  ir(J  • K<_Sh(K>£l) 


/*"  ^ Seo 

'L  / dK  Jx(Kp)  Ep(s,K,d)  — 


(H-5) 


(H-6) 


IK  J^Kp)  eQ  L~1(Ep(s,K,s)  J-) 
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/t~t' 


dt'  / dK  KJ 


/"  “ “i(^)  V**>  Ji«» 


(va) 


V? 


+ Hsdav(t,0,d) 


(H-7) 


from  which 


-I. 


®$,dev^t,p,<^  % yady 


C dt" 

t— t1* 

/*  dt'  f(t-t*-t") 

dt"  9a(c 

Miw*2j 

4 j 

6 

\ 

pod  / 

U£ 

It 


_[ya(p2+c2t’2)3/(-4t”)  y ,yccpt\ 

- e I1<~TPr~) 


if  a <<  max 


l Vo/ 


(H-8) 


-I  /-t 

«agr  A dt"/ 


dt*  f(t-t,-cM) 


d 9 /0 i i»t" 
dip  °l~2 
\ pod  , 


it1  / yo  \ 3 (-yoa^)/ (4t")  /At"  _ 2\ 

4/y  E \ 2t")  \po  3 / 


o“o 


if  a >>  max 


(H-9) 


Thus,  the  5-respoase  for  a * a 3(t  ~ “)/p  is 


ix  J<,(Ka>  Ji«5)  e 


(-ITe M)/(„c) 


(H-10) 


it  becomes 


d . /«,!«"■ 


JuaCp^+c  (t-t")  ; ]/ (~4t")  I^jpoc0(t-t>,)j 

/n.  iTTt  \ 1 via  Viocp  f 

2t  7 


If  *.  *.  »*  > 


^^(-MOC2!^)/ (4t) 
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jo  2a  ! 

fd  A /< 

t>|  i^t  \’ 

( uc  od  I 

[df  64 | 

J|  ] 

pod  /- 

.{-pop  )/(4t)  pgp 
: 4t 


(B-ll) 


Here  a condition  on  the  time 


e(iiac  t)/  4 ^ x 


(H-12) 


has  been  used  to  validate  the  third  equality  in  (H-ll).  For  cases  of 

practical  interest,  (H-12)  is  amply  satisfied  for  times  up  to  the  peak 

-1  2 

diffusion  time  t 10  pad  at  z * d,  because  usually  we  have 

2 

(peed)2  »!•*-♦•  1.4  x 10i-  x / u,  x — 2_  x -JL_\  >>  i 

' 2x10'  10“-*/ 

(H-13) 

(5) 

On  the  other  hand,  if  a » p the  , becomes 

$ *dev 


H 


(«) 

$ ,dev 


(t,p,d) 


dt" 


_L  r-pq  uo(a2+c2(t-t")2)/ (-4t") 

at  j_2t" 


T / pqac(t  - 
o\  2t" 


„ ^ to  2a  T d „ \1  . (-yoa")/ (4t)  yap 

4'  ~ 


uoa 

4t 


(H-14) 
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